Quantum fluctuations in the shape of exotic nuclei 



o 
o 

Oh' 

<D 



o 

=5 
S3 



o 
o 



X 



J. M. Yao 1 ' 2 - 4 , J. Meng 2 - 3 , P. Ring 4 , Z. P. Li 2 , and K. Hagino 5 
1 School of Physical Science and Technology, Southwest University, Chongqing 400715, China 
2 State Key Laboratory for Nuclear Physics and Technology, 
School of Physics, Peking University, Beijing 100871, China 
3 School of Physics and Nuclear Energy Engineering, Beihang University, Beijing 100191, China 
4 Physik-Department der Technischen Universitat Miinchen, D-85748 Garching, Germany and 
° Department of Physics, Tohoku University, Sendai 980-8578, Japan 
(Dated: September 9, 2009) 

Quantum fluctuations concerning the shape of nuclei are treated within the framework of covari- 
ant density functional theory. Long range correlations beyond mean field are taken into account 
by configuration mixing of wave functions with triaxial shapes and the restoration of spontaneously 
broken rotational symmetries through three-dimensional angular momentum projection. The con- 
troversial nucleus 16 C is treated as an example and it is found that its ground state has a triaxial 
shape but with large shape fluctuations. They are of crucial importance for a proper description of 
the spectroscopic properties of such nuclei. 
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Atomic nuclei are highly correlated and strongly inter- 
acting many-body systems of mesoscopic character. On 
the one side they are so small that quantum effects are 
dominant, on the other the number of their constituents 
is finite and large enough so that classical concepts such 
as their shape and orientation play an important role for 
an understanding of their structure. Only in heavy nuclei 
are the deformation parameters well defined quantities. 
In transitional nuclei, and in light systems, quantum fluc- 
tuations have to be taken into account. 

The rapid development in radioactive nuclear beam fa- 
cilities and gamma ray detectors have in recent years al- 
lowed one to study exotic nuclei far from stability and 
many new phenomena have been observed and predicted 
in this context. In order to investigate the shape de- 
grees of freedom it is not sufficient to only study the 
ground state properties of these nuclei, one also needs 
information about the spectroscopy of the low-lying ex- 
cited states. Therefore, presently much interest is focused 
on the measurement of the energies of the first 2 + , or 4 + 
states and of the reduced transition probabilities (B(E2)- 
values) from the first 2 + (2^~) to the ground state (0^~) [l|. 
These are fundamental quantities which reveal rich infor- 
mation about nuclear shapes and shell structure. 

Recently, the structure of the nucleus 16 C having two 
neutron particles and two proton holes above the dou- 
bly magic nucleus 16 O has become a very interesting and 
challenging topic. In lifetime measurements using the re- 
coil shadow method (RSM) |2j an anomalously hindered 
B(E2 : 2+ -> 0+) value of 2.6(9) e 2 fm 4 was found, con- 
tradicting the low excitation energy of E(2+) = 1.766 
MeV (see Ref. Q for an earlier measurement). Further- 
more, inelastic proton scattering 0, HI indicated a large 
quadrupole deformation (3 PP ~ 0.47(5). The quenched 
B(E2) value, combined with the large nuclear deforma- 
tion, led to the suggestion that the neutron motionplays 
a predominant role in the 2f state of 16 C 0, S 0|, al- 



though the structure of 16 C has not yet been fully un- 
derstood. The more recent lifetime measurement of the 
2^ state using the recoil distance method (RDM) after 
a fusion-evaporation reaction gave a similar, but slightly 
larger, B(E2) value of 4.15(73) e 2 fm 4 0. 

On the theoretical side, very different models have been 
used to describe the quenched B(E2) value in 16 C, in- 
cluding antisymmetrized molecular dynamics ( AMD) [8[ , 
three-body models 0, [l(J EH and the shell model . 
In AMD, the unusually small B(E2) value derived from 
the lifetime measurement with RSM was interpreted as 
the coexistence of an oblate proton and a prolate neutron 
shape. In the three-body models or in the shell model, a 
careful adjustment of an effective charge or modifications 
of the Hamiltonian were required to reproduce the data. 

The self-consistent mean field approach derived from 
a global density functional theory (DFT) provides a 
vivid way to study macroscopically defined quantities, 
e.g. ground state energy, nuclear radius and deforma- 
tion. Exotic shapes of 16 C have been investigated using 
non-relativistic and relativistic DFT calculations, con- 
strained by deformation parameters [3, 14 1. However, a 
very flat energy surface has been found both on the pro- 
late and oblate side indicating that this is a transitional 
nucleus. Therefore, substantial effects from configuration 
mixing connected with shape fluctuations are expected 
which play an important role for quantities of a quan- 
tum nature, e.g. , discrete energy spectra and transitions 
probabilities [15j | . In order to understand such matrix el- 
ements, e.g. B(E2) values, a microscopic approach going 
beyond the mean-field level is required which is able to 
treat shape quantum fluctuations properly. 

In the framework of DFT, shape fluctuations and angu- 
lar momentum projection have been treated using non- 
relativistic density functionals of Skyrme \m hj and 
Gogny [18], as well as covariant functionals [13]. In all 
of these cases, however, intrinsic triaxiality was neglected 
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which will be found to be crucial for a reproduction of the 
data. Only very recently, a mixing of all five quadrupole 
degrees of freedom has been attempted within the con- 
text of triaxial Skyrme calculations [2(| • 

The present investigation is based on a similar idea us- 
ing covariant density functional theory (CDFT). Based 
on Lorentz invariance, this method connects the spin and 
spatial degrees of freedom of the nucleus in a consistent 
way . Numerous investigations have shown that the ex- 
perimental data for the ground and excited states can be 
nicely interpreted within a relativistic framework. Re- 
cently, we have developed three-dimensional angular mo- 
mentum projection (3DAMP) for relativistic point cou- 
pling models to incorporate correlations related to the 
restoration of broken rotational symmetries 2 if . This 



concept has now been further extended to include fluc- 
tuations for triaxial deformations within the framework 
of the Generator Coordinate Method (GCM). 

In this letter, we use this method to study the effects of 
quantum fluctuation for triaxial shapes in the controver- 
sial nucleus 16 C. It turns out that the ground state of 16 C 
corresponds to a triaxial shape having large shape fluc- 
tuations along the 7 degree of freedom, which describes 
the triaxiality. Such a new picture for the shape of 16 C 
is essential to reproduce quantitatively the experimental 
B(E2 : 2+ -> Of ) value of Ref. that is consistent with 
the results from inelastic {p,p') scattering 

The set of deformed intrinsic wave functions \q) with 
the quadrupole deformations in the Hill- Wheeler co- 
ordinates q = (/3,7) [IH is generated by constrained 
RMF+BCS calculations using the parameter set PC-F1 
of Ref. [22j ■ Further details are given in Ref. 21 1 . In par- 
ticular, the strength parameters of the zero-range pairing 
force are V n = 308 and V p = 321 MeV-fm 3 . 

In the left panel of Fig. Q] we plot the potential energy 
surface (PES) in (3--f plane obtained with these mean 
field calculations. This PES is rather soft, which is con- 
sistent with the results found in Ref. Of course, for 
all the non-spherical points with (/3 ^ 0) the intrinsic 
wave functions |/3 7^ 0,7) have a certain orientation and 
rotational symmetry is spontaneously broken. To restore 
this symmetry, 3DAMP [2l[ is carried out and the right 
panel of Fig. [T] shows the projected PES for J T = + . 

Of particular interest is the occurrence of an obvious 
triaxial minimum with (3 = 0.6, 7 ~ 20°. To demonstrate 
the shapes at the various points in the (/3-7)-plane, in 
Fig. [2] we plot the corresponding intrinsic density distri- 
butions for neutrons and protons in the y-z plane inte- 
grated over the x-coordinate. Fig. [2^ refers to the axially 
symmetric minimum in the unprojected PES of Fig. [Ij, 
and Fig. [2)3 corresponds to the minimum in the projected 
PES of Fig. QJ). Here we observe rather different defor- 
mation parameters for neutrons (/3„ = 0.69, 7„ = 15.48°) 
and protons (/3 P = 0.46, 7 P = 31.47°), which shows us an- 
other kind of decoupled structure for the density distri- 
bution of neutron and proton, different to that indicated 
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FIG. 1: (Color online) Potential energy surfaces in /3-j plane 
for the nucleus 16 C obtained by triaxial relativistic mean-field 
calculations (left panel) and with projection onto angular mo- 
mentum J = after the variation (right panel). The energy 
gap between two neighbor contour lines is 0.5 MeV. 



in Ref. [8]. A careful analysis shows that the difference 
between neutron and proton deformation is due to the 
special 2p-2h configurations v(ld 5 /2) 2 <8> ^(^3/2) 2 ■ For 
large prolate deformations ((3 > 1.2), the deformation 
driving orbit ld 5 / 2 dives into the Fermi sea of protons as 
well. In this case, neutrons and protons have almost the 
same deformation and the decoupled structure ceases. 

Although there is an obvious minimum in the projected 
PES shown in Fig. [TJ this minimum is still very soft along 
the 7 direction. One has, therefore, to allow for superpo- 
sitions of wave functions with different deformations and 
to perform a GCM calculation in the full f3-j plane (for 
details see Ref. (Ell) 
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including correlations due to restoration of broken sym- 
metries and fluctuations of the deformation coordinates. 

To illustrate the importance of shape fluctuations in 
the 7 degree of freedom, in a first step, we restrict 
ourselves to axially symmetric shapes and choose the 
generator coordinates as /3=0, 0.1, • • • , 1.5 for 7=0°, and 
7=180°. As the result, the full 3DAMP+GCM calcula- 
tion is simplified into a 1DAMP+GCM calculation and 
all components with K ^ vanish. In Fig. [3^ we plot 
energies and average quadrupole moments of the lowest 
GCM states for each angular momentum + , 2 + , 4 + to- 
gether with the mean-field PES and the projected PES 
as a function of (3. The corresponding probability dis- 
tributions \g'a{(3)\ 2 are shown in Fig. [3b where <?^(/3) is 
the solution of the Hill- Wheeler- Griffin equation in the 
"natural basis" (for details see Ref. [15j]). 

Fig. [3] shows two minima on the projected PES with 
J" = + , one on the prolate and one on the oblate side 
having similar probabilities. Therefore, 1DAMP+GCM 
calculations lead to a nearly vanishing quadrupole mo- 
ment for the ground state with a rather small B(E2 : 
2 t Q t) = 0.11e 2 fm 4 . According to Fig. [TJ it is es- 
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FIG. 2: (Color online) Intrinsic density distributions of protons and neutrons in the y-z plane (the z-axis has been integrated 
over) (a) for the minimum of the mean-field potential energy surface (PES) of Fig. []J, at j3 = 0.4, 7 = 0°; (b) for the minimum 
of the projected PES (0+) in Fig. [lb with (3 = 0.6,7 = 20°; (c) for the deformation of the projected GCM state 0^ with 
(p) =0.44, (7) =24.27°. 



scntial to include triaxial states and to carry out a full 
3DAMP+GCM calculation. The probabilities in Fig. [3b 
decay quickly for (3 > 1.1. It is therefore justified to 
choose triaxial states with < f3 < 1.1 and < 7 < 60° 
having steps A/3 = 0.1 and A7 = 10° respectively. 
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FIG. 3: (Color online) Results of 1DAMP+GCM calcula- 
tion with axially symmetric intrinsic states: (a) energies and 
average quadrupole moments of the lowest GCM states for 
r = + , 2+, 4+ in 16 C, together with the mean-field (dotted 
hne) and projected energy curves, (b) squares of collective 
wave functions 
states. 



|<?a(/3)| for the corresponding lowest GCM 



Fig. [4] shows contour plots of the probability distribu- 
tions in the plane resulting from 3DAMP+GCM cal- 
culations for both the ground state (0]*~) and the first ex- 



TABLE I: Average deformation parameters (/3), (7) for the 
ground state (0^) obtained from 3DAMP+GCM calculations. 
The quantities /3mhi° 1 7min an< A /^mini 7min are the deformations 
at the minima of the projected (J = 0) and mean-field poten- 
tial energy surfaces. 





(P) 


(7) 


/3 J=0 

Mm 111 


"V J=0 
7min 


oMF 
rmin 


MF 
7min 


neut. 


0.50 


21.41° 


0.69 


15.48° 


0.48 


0° 


prot. 


0.34 


31.30° 


0.46 


31.47° 


0.27 


0° 


total 


0.44 


24.27° 


0.60 


20.00° 


0.40 


0° 



cited state (2^~). The ground state has a triaxial structure 
with almost uniform probability along the 7-direction 
which indicates an obvious quantum shape fluctuation in 
7 direction. For J = 2 we have if-mixing and we there- 
fore show the two distributions for K = and K = 2 
separately. Both are concentrated along the axially sym- 
metric configurations with K = on the prolate and 
K = 2 on the oblate side. However, the K = part ex- 
hausts 93.1% of the norm so that the nucleus has a strong 
prolate deformation in the (2^~) state with (3 ~ 0.6. 

For comparison with the mean-field results, we define 
the quadrupole moments for the ground state 

(Q 2M ) = I d 2 q\glz°M\ 2 q2Mi M = 0,2. (2) 



The deformations (j3) = 47r v /(Q 20 ) 2 + 2(Q 22 } 2 /(3iVi?g) 
and tan (7) = V2(Q 2 2) I '(Q20) shown in Tab. U result 
from 3DAMP+GCM calculations. We observe a decou- 
pled structure, where the average neutron deformation 
(/?)„ = 0.50 is obviously larger than that of protons 
((3)p — 0.34. This can also be seen in Fig. [5] where 
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f (r) = J d 2 q\gizl{q)\ 2 Pr{r ]q ) (3) 
are the density distributions of neutrons and protons for 
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FIG. 4: (Color online) Contour plots of the probability dis- 
tributions in f3 and 7 from 3DAMP+GCM calculation for the 
ground state and the first excited state in 16 C. 
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FIG. 5: (Color online) The lowest energy levels of angular 
momentum J = 0, 2, 4 in 16 C. The B(E2) values are in units 
of e 2 fm 4 . The experimental data are taken from Refs. @,@]- 



the projected GCM state (0^). 

To demonstrate the importance of shape fluctuations 
in 7 for spectroscopic properties, we plot in Fig.[5]thc low- 
est energy levels having angular momentum J = 0, 2, 4 in 
16 C obtained from the ID and 3D AMP+GCM calcula- 
tions and compare them with experiment The inclu- 
sion of triaxial states does not change the energy spec- 
trum too much, but it improves the electric quadrupole 
transitions significantly. Only with a proper treatment 
of the shape fluctuations in 7, are the calculated B(E2 : 
2f — > 0i) = 6.50 c 2 fm 4 in agreement with recent data. 

The predicted ratio -R4/2 = E4/E2 — 2.63 is a little 
larger than the experimental value R4/2 — 2.35. Both are 
close to the R4/2 = 2.50, a typical value for 7-soft nuclei. 
The predicted deviations of R4/2 from experiment go in 
the direction of a rotor (R4/2 = 3). This has its origin in 
the increasing prolate deformation as angular momentum 
increases to J = 2 and J = 4. As observed in many 
calculations using projection after variation (l6l . 17. 18. 



[l9j the calculated spectrum is systematically stretched. 
It would be far more compressed by the inclusion of a 
cranking term in the mean- field calculation [23| . 

In summary, starting from covariant density functional 
theory and including additional correlations by restor- 
ing symmetries and by configuration mixing of triaxial 
shapes we have studied the influence of quantum shape 
fluctuations on spectroscopic properties in the exotic nu- 
cleus 16 C. In contrary to earlier investigations using other 
methods we find that the ground state has a triaxial 
shape with large shape fluctuations in 7. Such a novel 
picture for the shape of 16 C is essential to reproduce 
quantitatively B(E2 !) values derived from a recent life- 
time measurement based on the RDM method. 
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